Persistent currents in mesoscopic rings with a quantum dot 
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Using the Anderson model in the Kondo regime, we calculate the persistent current j in a ring 
with an embedded quantum dot (QD) as a function of the Aharonov-Bohm flux $ for different ring 
length L, temperature T and broadening of the conduction states 5. For T — S = and L S> f, 
where £ is the Kondo screening length, Lj tends to the value for a non interacting ideal ring, while it 
is suppressed for a side coupled QD. For any L/£, Lj is also suppressed when either T or 8 increase 
above a fraction of the level spacing which depends on $. 



Electron transport through a quantum dot (QD) has 
been a subject of great interest in the past few years. The 
progress in nanofabrication made it possible to use QD's 
as ideal realizations of the Kondo effect, which is one 
of the most exciting and studied problems in condensed 
matter physics, ft consists in the screening of an impurity 
spin by a cloud of conduction electrons of radius £ ~ 
Tivf/Tk and energies ~ 2Tk around the Fermi energy ef, 
where Tr- is the Kondo temperature and vf is the Fermi 
velocity. The energy scale Tk is experimentally accessible 
in different ways, like the width of the peak in linear 
response conductance through a QD and its temperature 
dependence Q] . Instead, a direct measurement of £ does 
not exist so far. 

Several interesting experiments were performed re- 
cently in an Aharonov-Bohm geometry, in which a ring 
containing a QD is threaded by a mag netic flux ||-||. 
Phase coherence along the ring has been demonstrated. 
The persistent current j in these rings, and in rings side- 
coupled to a QD has been studied theoretically [p|-pd|. 
However, basic results of these works contradict each 
other and an accurate method to calculate j for any 
L/£ has not been developed. Exact Bethe ansatz results 
were known for L 3> £ and chiral electrons 0. These 
were extended to electrons moving in both directions || . 
The precise geometry of these calculations was explained 
recently Perturbative renormalization group (RG) 

calculations have established the form of for L <C £ 
and L — ► oo @]. A change in the dependence of j with 
magnetic flux is expected between the regimes L -C £ and 
L ^ £ dhll- Thus, measurements of j would provide a 
way of detecting £ 0. Since the average level spacing 
is D = 2ttTivf/L, the condition L ~ £ is equivalent to 
D/2ir~T K . 

Our purpose is to describe j($>) accurately through the 
crossover region, and to consider the effects of tempera- 
ture T and finite level width 6 of the conduction states 
for the first time. We relate j with the one-particle Green 
function at the dot and the latter is calculated using an 
interpolative perturbative approach (IPA) |l3|-|l5|]. For 
small L and T = 6 = 0, we also calculate j using numer- 
ical exact diagonalization (ED) finding excellent agree- 
ment with the IPA. 



The Hamiltonian for the embedded dot is: 

L-l 

II X (* 3 -e^4 +1<r c,v+H.c. 

>r,3=0 

+E d ^n da + U n d ^n di , 



(1) 



with clct = coer and n da = c^coa- The phases cpj depend 
on the choice of gauge and satisfy ^ ■ (fj — <f>, where 
$/ic/27re is the magnetic flux threading the ring. The 
hoppings are all equal except those involving the QD: 
tj = tR (right) for j — 0, tj — tr, (left) for j = L — 1, and 
tj = t otherwise. The total current flowing between sites 
I and I + 1 is: 



e d{H) 
h dipi 



le \ ^ , 



Cl a - H.C. 



(2) 



The expectation value entering Eq. (Q) is given in terms 
of Green functions ((ci a ; cj CT ))- These in turn can be ex- 
pressed in terms of the QD Green function G da (ui) = 
((coo-; Coa))^ using equations of motion. Choosing <pj = 
for < j < L — 1 and ipa = fL—i = 7r/2, we obtain, after 
some algebra for ^ I ^ L — 1: 
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tt R t L sin $ ^(-I) n ^ p(/, n)p(l + 1, n) 



/ du>f(u))lm 



G da (uj + irj) 



(w + irj— e n )(uj + irj - e n <) 



Ttnl 7rn Ttn 

p(l, n) = sin — — sin — ; e n — —2t cos — id. (3) 

Lj L L 

Here ~^2 n runs over all integers 1 < n < L— 1, while Yl' n ' 
is restricted to integers of opposite parity to that of n, 
f(to) is the Fermi function, rj is a positive infinitesimal, 
and e n are the eigenergies in absence of the QD, allowing 
for a finite broadening 5 [^2|. Eq. (j|) is exact. However, 
we use an approximate G da - Within the precision of the 
numerical integration, we have verified that the resulting 
ji is independent of I, as it should be because of conser- 
vation of the current. Therefore we drop the subscript in 
what follows. 
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We calculate Gda using a self-consistent IPA based on 
perturbation theory in U up to second order , gen- 

eralized to allow spin dependence: 

where G^ a is the Green function for U = and Ed re- 
placed by an effective energy e" ^ determined selfconsis- 
tently 

HMY 1 -ES 

?Jcr = (nda), and: 

s / s = M 1 - n^)sj 2) 

n°(l - n°) - [(1 - ^)C7 + £ d - E^tf-^ ' 

where n° is the expectation value of n d(T calculated with 
G da and is the ordinary second order correction to 
the self energy, calculated from a Feynmann diagram in- 
volving the analytical extension of G da (ui) to Matsubara 
frequencies [|l2| : 

^(m,T) = U 2 Tj2G° d(T (i"i - w m ) x {w m )\ 

m 
a 

The resulting Gda is not only valid up to U 2 , but it is also 
exact for a decoupled dot (£t = tn = 0), and reproduces 
the leading term for lu — > oo p^ ]. We determine by 
imposing that n£ = n<7 for both spins. 




FIG. 1. Current in units of et/% as a function of magnetic 
flux for U = 2, tn = 0.4 and different values of L. Unless 
otherwise indicated N = L. 



We take t = 1 as the unit of energy and keep = 
— U/2 and it = i/j. This allows us to exploit electron- 
hole and reflection symmetry in some cases. As a ba- 
sis for our study we choose U = 2 and = 0.4. For 
L — ► oo and £f = 0, this leads to a resonant level width 
A = 0.32 p|| (neglecting its energy dependence). The 
ratio U / A = 6.25 is large enough for the system to be in 
the Kondo regime of the model, but low enough to en- 
sure the validity of the IPA p5[ . The impurity spectral 
density Pda(u) shows three peaks at Ed, Ed + U and £f 
characteristic of the Kondo regime. For 6 = and finite 
L, pdcr(w) consists in a set of delta functions. Using [ p^[ : 

the above parameters lead to Tk ~ 0.05 and £ ~ 40. We 
begin showing the results for T = 5 = 0. They depend 
drastically on the parity of the number of particles N. 
For even N, the results for odd N/2 are essentially the 
same as those for N ± 2 shifting $ by 7r. Then we re- 
strict to cither odd N or A/4 integer. Using reflection 
symmetry around the QD (cj a — ► C£_j CT ), one realizes 
that j($) = — j(— Thus, it is sufficient to represent 
j in the interval < 3> < ir. Fig. 1 displays the evolu- 
tion of as a function of ring size for even N. For 
L = 8 and A = 7, 8 we also compare the IPA results with 
those obtained using j($) = -(eL/K)dE(3>)/d$, with 
the ground state energy E calculated by ED. The maxi- 
mum deviation between both results takes place around 
$ = 0.2 and is below 0.05et/h. For L — > oo and A even, 
i($) converges to the non-interacting ideal result, as pre- 
dicted by an analysis of the strong-coupling fixed point 
of RG @ . This supports the validity of the IPA results 
for all L. Our results also agree qualitatively with those 
of RG for L <C £, but disagree with those of Ref. ||. 
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FIG. 2. Current as a function of L for two sets of parame- 
ters. 

To study the scaling properties in the dependence with 
L, we have calculated j{^/2) for 8 < L < 1000 and 
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two sets of parameters: U — 2, = 0.4 as before, and 
U = 0.5, <l = 0.2. The new choice should reduce A by 
a factor ~ 4 [[l5| , keeping then nearly the same value of 
U/A, and resulting in ~ 4 times smaller Tk (see Eq. (El)), 
and <~ 4 times larger £. In fact, as shown in Fig. 2, scaling 
L by a factor 3.79, both functions Lj(Tr/2,L) practically 
coincide, confirming that Lj is a universal function of 
L/£, as expected for L ^> 1 in the Kondo regime [Q. For 
L ~ 8 and f7 = 2, is shifted a little bit upwards from 
the universal curve. For U = 0.5, using £ = 151.6 as 
derived from Eq. (Q) , the results for 8 < L < 240 fit with 
negligible errors on the curve: 

Lj(n/2)h/te = 1.53 + 0.241n(L/£). (5) 



screening. Also, j (<&) displays positive and negative val- 
ues in the interval < $ < it, suggesting a tendency to 
periodicity in ir instead of 2ir. This periodicity is exact for 
N = L odd, Ed = —U/2 and tt — t^: the electron-hole 
transformation Cj a — > Cj CT , ^ — > y>j + 7T, maps H{— <f>) 
onto if(<f> + 7r) and j(-$) onto -j(<f> + 7r) (see Eqs. (|l|) 
and (||)). Combining this with = — j(— <f>), one has 

= + 7r )> and J^ 71 "/ 2 ) = 0- In Fig- 3 we snow 
j($) in two of these cases with U reduced to 1.5 to be 
able to obtain self consistent solutions up to L = 37. 
We estimate £ ~ 25. According to RG arguments, as 
L increases, tends to the result of a fictitious non- 
interacting system with spin dependent Ep to allow for 
odd N. Our results are consistent with this. s z decreases 
from ~ 0.37 to ~ 0.15 as L increase from 9 to 37. 
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FIG. 3. Same as Fig. 1 for odd N and U = 1.5. 

Using ED, we find that in general for given sf, the 
states with odd N are less stable than those of even N. 
In particular, it is not possible to find a unique value of 
Ef for which some odd N are stable for all <F Even N 
with odd (even) N/2 are favored near $ = ($ = 7r), 
for which j{$>) has a small amplitude. For example, for 
L = 12 and ep = —0.7, N decreases from 10 for $ = to 
8 for <I> = 7r, jumping from 10 to 9 near $ = 0.167T, and 
from 9 to 8 near $ = 0.827T. As a consequence j($) is 
discontinuous and small in magnitude (|i(3?)| < 0.5et/K). 
The IPA for odd N was applied imposing a given N (al- 
though it might correspond to a metastable state) and 
adjusting £j? as a function of A technical difficulty is 
that for large U and L we could not find the selfconsis- 
tent solution. This is related to the fact that for finite 
i, n° and n a are discontinuous functions of the £%tf and 
it is not always possible to satisfy = n a . As seen 
in Fig. 1, and the above mentioned case for L = 12, 
j($) is strongly suppressed for odd N and small L. This 
is related to a partial suppression of the Kondo effect. 
For L = 8, the expectation value of the impurity spin 
s z = (nj — Tii)/2 ~ 0.4, indicating only a small partial 




d>/jt 

FIG. 4. Same as Fig. 1 for different temperatures and (a) 
U = 2, t L = 0.4, L = 800, (b) U = 0.5, t L = 0.2, L = 8. 

We have calculated the T dependence of j($) for even 
and several In Fig. 4, we show two cases: (a) 

T K ~ 0.05, L/£ - 20 and (b) Tr- - 0.012, L/£ - 1/20. 
From the known results for the conductance through a 
QD jl],[l5| one would expect Tk to be the relevant scale 
for the T dependence. However, it is a fraction of the level 
spacing D = 2ttTlvf/L in both cases. This is easy to un- 
derstand in case (a): for T <C Tk and L ^> £, the physics 
is still dominated by the strong coupling fixed point of 
RG for which the model reduces to a non-interacting one 
0. In turn, the conductance of the latter is strongly 
reduced for T ~ D <C Trt. In contrast to the infinite sys- 
tem, for which the Kondo peak in pd^UJ) looses approxi- 
mately half its intensity for T ~ Tk jl5) , when T K < -D 
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(L <C £), Gda is practically not modified with increas- 
ing T, until T reaches a sizeable fraction of D. As a 
consequence for all L the scale for the T dependence is 
~ D/b, but depends on $. Note that these arguments 
are independent of the parity of N. The last argument 
is in agreement with recent RG calculations, which show 
that for even N and Tk <C D, the screening of the local- 
ized spin takes place at T ~ D, while as discussed above, 
the screening is only partial for odd N [Q. In other 
words decreasing L, the Kondo effect is enhanced for N 
even and inhibited for N odd. 




FIG. 5. Same as Fig. 4 (b) for T — and several values of 



6. 



So far, we have considered S = and G da (uj). However, 
real mesoscopic systems are coupled to a charge reservoir 
which determines ep and discrete electronic states be- 
come resonances. The effects of this coupling have been 
studied recently |I^,[l7]]. The coupling to the conduction 
eigenstates (or to the impurity) is essential to explain the 
line shape in the projection of the Kondo effect to a re- 
mote location (mirage effect) jL2j. In the case D <C Tk, 
it is clear that j is strongly reduced when 8 reaches D. In 
the opposite case Tk -C D, for N even (corresponding to 
some e n near ep), p da (oj) near ep evolves with increas- 
ing 8 from two delta functions at both sides of ep to one 
peak centered at ep [fL2|. The corresponding change in 
i($) is shown in Fig. 5. The effect of the broadening of 
the levels is similar to that of increasing T, but a little 
bit weaker and more evenly distributed in $. 

We discuss briefly the side dot. The Hamiltonian is: 

L-l 

H = - £ (^4+1^° + H - c ') - td^2(4^ + H.c.) 



",3=0 



+Ed^n da + Un d ^i di , 



(6) 



where now ip = $/£ and n da = $ a d a . Proceeding as 
before, for 5 = the current becomes: 



2te 



{-27rL^sin(fe + ¥ p)/( 



£k) 



k 

+ < 'EE sin ^ +A; '^ + 1 )-^ 



k W 



/ duf(u>) Im 



(w + ii] - e k ) (w + if) - e k > ) 



}; 



-2t cos(fc + ip), 



(7) 



independently of I . The expressions for G d a are the same 
as before with V n replaced by V k = 1/VL. In the limit 
t d — > 0, L — > oo, the resulting resonant level width A is 
the same as that of an embedded dot with tz = in = td/2 
M|. Then, we take t d — 0.8, N = L + 1 even and other 
parameters as in Fig. 1 to study the dependence of j s (<&) 
on L. The result is shown in Fig. 6. In agreement with 
Ref. j s ($) is suppressed for large L. Note that the first 
term of Eq. (^) is the non-interacting one, which is of the 
order of 100 times for L = 799. The fact that the 

small j s (&) comes from a near cancellation of two terms, 
and only the second one is approximate, further supports 
the validity of the IPA. 




FIG. 

U = 2, 
L. 



6. Current for the side dot as a function of flux for 
td = 0.8, T — S = 0, N = L + 1 and several values of 



We have calculated the persistent current j( < I > ) in a 
ring with a QD using an interpolative perturbative ap- 
proach. The method is accurate enough to describe 
for all £/£■ The universal dependence of Lj with L/£ is 
displayed. The energy scale for the dependence on tem- 
perature and broadening of the levels is a fraction of the 
level spacing D instead of Tk- For the side dot j s {$) is 
small and decreases with L. For both systems and odd L, 
the period of j(<I>) is 7r in half filled symmetric rings. Be- 
cause of their larger stability and larger currents, states 
with even number of electrons N in embedded dots seem 
the most interesting experimentally. 
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